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ABSTRACT 



X 



Lyman-break galaxies are now regularly found in the high redshift Universe by searching for the break in the galaxy spectrum caused 
by the Lyman-limit redshifted into the optical or even near-IR. At lower redshift, this break is covered by the GALEX UV channels and 
small samples of z ~ 1 LBGs have been presented in the literature. Here we give results from fitting the spectral energy distributions 
of a small sub-set of low redshift LBGs and demonstrate the advantage of including photometric points derived from HST ACS slitless 
grism observations. The results show these galaxies to have very young, star forming populations, while still being massive and dusty. 
LBGs at low and high redshift show remarkable similarities in their properties, indicating that the LBG selection method picks similar 
galaxies throughout the Universe. 
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1. Introduction 

To understand galaxy evolution we need to understand the stel- 
lar properties of galaxies at high redshift. These galaxies include 
samples of mainly Lyman-break galaxies (LBGs; e.g. Steidel et 
al. 1996, 1999, Pettini et al. 2001, Shapley et al. 2003, Bunker et 
al. 2004, Ouchi et al. 2004, Burgarella et al. 2007), but also other 
types of galaxies such as Lya emitters (LAEs; e.g. M0ller & 
Warren 1993, Fynbo et al. 2002, Venemans et al. 2007, Nilsson 
et al. 2009a) or very red and massive galaxies (DRGs, EROs; 
Franx et al. 2003, Daddi et al. 2004). One common method 
of studying the stellar populations of these galaxies is by fit- 
ting their spectral energy distributions (SEDs) with spectral tem- 
plates, thereby determining luminosity-weighted properties such 
as stellar ages, dust content, metallicities, and stellar masses. 

For Lyman-break galaxies, the last decade has seen a large 
number of publications determining the stellar properties using 
SED fitting. Two of the first publications, Shapley et al. (2001) 
and Papovich et al. (2001) presented results for LBGs at z — 
2 — 3. They found that LBGs at these redshifts tended to be 
young, with ages of a few hundred Myrs respectively, have 
masses around 10 10 M Q and significant extinction of the order 
of Ay = 0.7 — 1.2 mag. The star formation rates (SFRs) in their 
samples were a few times 10 M© yr _1 . SED fitting results at 
higher redshifts seem to indicate some evolution in these prop- 
erties. Verma et al. (2007) and Yabe et al. (2009) fit LBGs at 
z ~ 5.5 and find that they are typically younger (a few times 
10 Myrs), less massive (a few times 10 9 M Q ), and less dusty 
(Ay = 0.3 — 0.7 mag.). They also find that the higher redshift 
LBGs have higher SFRs, of the order 50 - 100 M yr~\ in- 
dicating that these are small galaxies in the process of building 
up their stellar mass. The evolution seen from smaller, younger 
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galaxies at z ~ 5 to larger and less vigorously star-forming 
galaxies at z ~ 2.5 does not seem to continue to even lower red- 
shift. Overzier et al. (2009) presented a sample of local "LBG 
analogues" based on UV bright galaxies in SDSS. These galax- 
ies have masses of 6 x 10 9 M Q and SFRs of ~ 10 M yr _1 , 
i.e. they are of similar mass and have similar SFRs to LBGs at 
z ~ 2.5. 

Whereas SED fitting remains a popular method to determine 
stellar properties of LBGs, potential caveats with this method are 
the lack of a good wavelength spread in, or numbers of, spectral 
points and low signal-to-noise in the observations. In this pa- 
per we aim to avoid both these potential sources of uncertainty 
by observing a sample of pre-selected, lower redshift LBGs, at 
z ~ 1 (Burgarella et al. 2007). We also focus especially on LBGs 
with associated HST/ ACS G800L grism spectra, since the spec- 
tra will, at these redshifts, completely encompass the Balmer 
and 4000 A breaks, thus allowing a very good sampling of these 
breaks. These breaks, in turn, are very effective indicators of the 
age of the stellar population and dust content (Hamilton 1985, 
Balogh et al. 1999, Kriek et al. 2006). In Kriek et al. (2006) it is 
shown that not only the redshift, but also the stellar properties are 
better constrained when using spectroscopic observations cover- 
ing the Balmer and 4000 A break of the galaxy, rather than just 
broad-band photometric points. Kauffmann et al. (2003) show 
that stellar properties can be determined solely on the basis of 
the 4000 A break strength and the Balmer absorption-line index 
US a- Slitless grism observations with the HST G800L grism are 
particularly well suited to study these breaks in high redshift 
galaxies; with low wavelength resolution, but high sensitivity 
and low background. Here, we present the results of fitting 15 
LBGs with grism spectra, and show that the stellar properties 
are very well constrained by these observations. 
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In Sect. 12 we describe the data used in the SED fitting, and 
in Sect.[3]the fitting method is described. The results are given in 
Sect.|4]and we end with a discussion in Sect. [5] Throughout this 
paper, we assume a cosmology with Hq = 72 km s _1 Mpc~\ 
f2 m = 0.3 and Qa = 0.7. Magnitudes are given in the AB sys- 
tem. 

2. Data 

2.1. LBG sample 

The sample of LBGs used here is drawn from Burgarella et 
al. (2007), who presented the selection of LBGs using GALEX 
observations in the redshift range of z — 0.9 — 1.3. Their sample 
includes 420 LBGs in the Chandra Deep Field South (CDFS). A 
sub-set of 32 of these lie in fields covered by public ACS slitless 
grism observations (Kummel et al., in prep.). The ACS obser- 
vations were made with grism G800L, with an effective wave- 
length range of - 5 500 - 10 500 A and a resolving power of 
66 at 6365 A for point-sources. All public ACS slitless grism 
observations have been reduced by the PHLAG pipeline in an 
effort to provide science-ready grism spectra of almost 32 000 
unique objects (Kummel et al., in prep., but see also Pirzkal et 
al. 2004, Xu et al. 2007, Straughn et al. 2009 for similar ob- 
servations). For this project, we specially extracted 1-D spec- 
tra from flux calibrated 2-D spectra for the 32 objects covered 
by grism observations, using a dedicated extraction scripQ. The 
extraction provides error and contamination estimates for each 
wavelength pixel. Of the 32 spectra, four were not possible to re- 
trieve due to their faintness or proximity to the edge of the CCD. 
Of the remaining 28, 15 were chosen for follow-up spectral fit- 
ting for having either a precisely known redshift from ground- 
based spectroscopy (13) or for having a very well-defined spec- 
tral break in the grism spectrum, allowing a redshift to be de- 
termined (2). The sample is defined in Table Q] Ground-based 
spectroscopic redshifts were collected from all available public 
surveys in GOODS-S. Object LBG_61554 had conflicting red- 
shifts from VVDS (z = 0.378) and the K20 survey (z = 1.087). 
In this case we retain the K20 redshift as we see a break in the 
grism spectrum which cannot be explained with the lower red- 
shift. 

No particularly strong selection bias is expected in this sam- 
ple, compared to higher redshift samples, since the original se- 
lection by Burgarella et al. (2007) is expected to select the same 
type of galaxies as at higher redshift, and since the overlap with 
the grism spectra is a purely random sub-selection. The only 
bias may occur in the final selection from 28 objects with grism 
spectra, to the 15 used in the final analysis. Here, objects with- 
out ground-based spectra, typically fainter than the ones with 
ground-based spectra, are rejected. However, this bias may be 
considered mild, or insignificant, since the range of UV lumi- 
nosities of the final sample of LBGs covers roughly the same 
absolute magnitude range, or even go a bit deeper, than the typi- 
cal range at high redshift. 

2.2. Broad-band photometry 

To fit the SEDs of these galaxies, we need to collect photomet- 
ric measurements from broad-band observations of the CDFS, 
and from the grism spectra. The photometry in the CDFS was 
done with aperture photometry on the images available in the 

1 See 

www .stecf.org/ archive /hi a/ i noli v_extr act ion . php 



GOODS-S HST (in filters F435W, F606W, F814W and F850LP, 
Giavalisco et al. 2004), the ESO ISAAC releases (Retzlaff et 
al. 2010), the ESO VIMOS £/-band (Nonino et al. 2009), and 
the Spitzer IRAC 3.6 /im band images. The HST images were 
convolved with a Gaussian function of FWHM 0.8" to match the 
ground-based data. The apertures in all images had a 2" diame- 
ter. One source (LBG.61983) was not covered by the GOODS-S 
imaging, and the photometry for this source (in B Viz) was found 
in the MUSYC catalogues (Gawiser et al. 2006). Photometry in 
longer wavelength filters than the 3.6 /im were not attempted as 
these correspond to > 1.5 fim in the restframe where the spectral 
fitting models of our code are not very reliable. 

In order to test the accuracy of the photometry, specifi- 
cally the agreement between photometry from different instru- 
ments, two tests were performed. Firstly, the broad-band pho- 
tometry was run through the Hyperz photometric redshift code 
(Bolzonella et al. 2000). Only the 14 galaxies with GOODS-S 
detections were tested, of which eight were assigned the cor- 
rect redshift within the la significance interval directly, for five 
the correct redshift was given as a secondary solution and only 
one (LBG_53929) was assigned an incorrect redshift (z ~ 2.3). 
Hence, the photometry does not seem to suffer from any large 
sources of systematic error. Secondly, five unblended stars were 
found in the images and the procedure to extract photometry was 
repeated. The photometry from these stars were fitted with the 
stellar library of Gunn & Stryker (1983). With the minimum 
error of 0.07 magnitudes used in the SED fitting here, the re- 
duced x 2 °f all fi ve stars were for the best fits ranging between 
1.79 — 3.10. This test proves that there are no negative effects of 
potential aperture corrections, and that the minimum error bars 
chosen are optimally set. 

2.3. Grism photometry 

To constrain the Balmer break in these galaxies, photometry was 
extracted from the ACS slitless grism spectra in the following 
way. Spectral flux points in slitless grism spectra are correlated 
on small scales, and the correlation length along the wavelength 
direction depends on the size of the object. This is also mani- 
fested in a change in resolution with the size of the object. For a 
point-source object, the resolution of the grism is 66 at 6365 A, 
and the resolution scales linearly with object size according to: 

where 0.11 is the intrinsic pseudo-slitwidth in arcseconds for a 
point-source, and PSW is the pseudo-slitwidth in the dispersion 
direction of the actual object (as given by the APERW ID keyword 
in the released spectra). The PSW is in turn related to the mor- 
phological parameters of the object (i.e. major and minor axis, 
and position angle compared to spectral extraction direction). 
From the resolution we can calculate the FWHM and a of the 
dispersion. Extracting photometric points at 3a intervals apart 
should then ensure virtually uncorrelated measurements of the 
SED of the object. Finally, for objects with multiple spectra at 
different HST/ ACS G800L roll-angles, photometric points were 
extracted with 3a intervals, but with slight offsets in wavelength, 
so that the extracted points sampled the entire spectrum more 
evenly. The extracted points were normalised to the broad-band 
photometry with the HST i or z magnitude, in order to eliminate 
effects caused by different aperture sizes, and then input into the 
SED fitting code. The magnitude used for the normalisation was 
excluded from the fit. Two examples of the SEDs are shown in 
Fig.H] 
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Table 1. LBGs used in the fitting. 
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IDs are from Burgarella et al. (2007); RA and Dec (in J2000) from the grism position; My is the isophotal magnitude from the GOODS-S 
catalogue v.2.0; FWHM are measured in the z band from the GOODS-S catalogue. The last columns give number of overlapping grism 
spectra on each source (with different orientation angles), the number of extracted SED points for each object, and an example ID for the 
grism spectra. Notes on redshifts give the source of the ground-based redshift: 1) VVDS (Le Fevre et al. 2005), 2) IMAGES (Ravikumar et 
al. 2007), 3) FORS (Vanzella et al. 2008), 4) K20 (Mignoli et al. 2005) and 5) VIMOS (Balestra et al. 2010). 
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Fig. 1. SEDs and best fit spectra of two galaxies, one with a good xf- an d one with a worse fit. Points are photometric points, from 
broad-band photometry (brown) and from the grism spectra (red). Error bars in the wavelength direction show the FWHM of the 
filter with which the point was observed. The grey spectrum is the best fit spectrum from the models, with green, yellow and blue 
lines showing the spectra convolved with the resolution profile for each galaxy for the 1 SSP, 2 SSP and constant SFR scenarios. 
LBG_46190(/e/f)hasaXr ~ 1 whereas LBG.53929 has a xl ~ 25 (.right). The brightest of the grism data points in the LBG.53929 
SED is due to a strong [Oil] emission line. 



3. SED fitting method 

For the fitting, the spectral models of GALAXEV were used 
(Bruzual & Chariot 2003). The actual fits were made with the 
Monte-Carlo-Markov-Chain (MCMC) code NisseFit (Nilsson et 
al. 2007, Nilsson et al. 2010). We have described this code in 
some detail in previous publications, and will here only give a 
brief summary. The code uses a Metropolis Hastings algorithm 
to produce a Markov-Chain in a multi-dimensional parameter 
space. At each iteration point, the "quality of fit" parameter (in 
our case x 2 ) 1S calculated by running the GALAXEV code with 



the given set of parameters and a Salpeter IMF. For single stellar 
population (SSP) models we add both nebular continuum and 
nebular emission lines based on age, metallicity and UV flux. 
These were calculated independently on a grid of the three pa- 
rameters using the SB99 (Leitherer et al. 1999) and MAPPINGS 
(Kewley et al., in prep.) codes. The output spectrum is also cor- 
rected for the effect of dust according to Calzetti et al. (2000) and 
redshifted. The model magnitudes in the different wavebands are 
extracted from the final spectrum and compared with the ob- 
served values to produce the "quality of fit" parameter, where 
we use a flat prior. 
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For the fits here, from restframe ~ 0.18 — 1.8 /im, three runs 
were made with constant star formation rate (SFR), one SSP or 
two SSPs. The reason to fit SSPs is that we are only able to add 
nebular lines to these models, and two SSPs allow for a more 
realistic scenario to be tested. Single SSPs are then also run in 
order to compare the results, as many other studies have been 
made with single SSPs. A test of the validity of the SSP star 
formation history is then to compare the results to those of the 
runs with constant SFR. In all cases the age, metallicity, dust Ay 
and stellar mass are free parameters. For constant SFRs, the SFR 
is also fitted, while for two SSPs we fit two ages and the mass 
fraction between the two populations. The minimum error on the 
broad-band measurements was set to 0.07 magnitudes to account 
for absolute calibration uncertainties. Redshifts were assumed 
to be those determined from the ground-based spectra. For the 
two without such values, the redshift was set to z = 1.0. Fits 
made with redshifts varied in increments of z — 0.02 between 
z = 0.9 - 1.1 all revealed similar results. Objects LBG_49138 
and LBG_53929 have very strong emission-lines in their spec- 
tra (see e.g. Fig. [1} that are fitted poorly by the SED code. In 
the first attempt to fit 2 SSPs, the respective reduced x 2 were 87 
and 24. Manually removing 2/1 points in the grism photometry 
improved the fits to \ 2 r ~ 8 for both galaxies. We report these 
results in the next section. 



4. Results 

4.1. SED fitting results 

NisseFit was run 10 000 times for each galaxy, and the first 2 000 
runs were excluded to reach the fits where the x 2 had converged 
to values near the minimum. The resulting runs then reflect the 
actual best fit distribution of the parameter probability space fit- 
ted. To study these probabilities, we collapse the runs into all 
possible 2-D combinations, and plot the significance contours. 
One example is shown in Fig. [3] It is clear that metallicity (in 
all fits) and SFR (in the constant SFR scenarios, not shown in 
Fig. cannot be determined accurately, and we ignore these re- 
sults hereafter. The parameters that are well determined are stel- 
lar mass, dust Ay and age, although the age fits in the 2 SSP 
fits suffer from a degeneracy at the older age when the young 
population age is very young. In other words, when the young 
population age is very young, the light of these stars will out- 
shine any older population, unless the mass fraction in the young 
population is very small (see Nilsson et al. 2010, for a longer dis- 
cussion on this). The ages and mass fractions for all the 2 SSP 
fits (the only runs with more than one age fitted) indicate that 
we are not able to identify an older population in any of the ob- 
jects, and only upper limits on the mass fraction in a potential 
old population may be inferred. 

To determine the best fit parameters in 1-D, a method similar 
to that of Nilsson et al. (2009b) was used, where each parame- 
ter is decided individually and independently from the other pa- 
rameters by finding the median of the distribution of runs and 
integrating to 68.3% of each wing to find the scatter. The best fit 
results for the runs with two SSPs are found in Table|2] Here we 
only give the results of the runs with two SSPs, as the results for 
the constant star formation rate and 1 SSP runs agree very well 
with those for two SSPs. We see that most objects are fitted well, 
with reduced x 2 S 3- Objects with worse x 2 can be explained 
by features in the SEDs, from emission lines or too small error 
bars on the photometry due to systematic errors. 



Table 2. Results from the SED fits for the LBGs. 
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Stellar properties from fitting 2 SSPs to the data. Columns give the 
age of the younger SSP, dust Ay, mass in stars, the minimum re- 
duced x 2 an d tne dust corrected SFRs as calculated from the UV 
flux of all the fits. We here only give the age of the younger SSP, 
as any older stellar population is totally obscured by this younger 
population. LBGs 49138 and 53929 have been fitted after removing 
obvious strong emission lines. The bad x 2 of 62522 suffers from 
odd features in the NIR ISAAC photometry. Both object 53929 and 
66235 lie on the edge of the HST imaging field and may thus also 
have larger systematic uncertainties on the photometry. Please note 
that the error bars in column six are the pure measurement errors 
on the U band photometry, and do not incorporate uncertainties in 
the extrapolation to 1500 A, or in the SFR conversion factor. 



4.2. Star formation rates 

From the U band photometry, probing restframe ~ 1600 — 
2150 A at z = 1, we can calculate the UV-derived star formation 
rates of the galaxies. Using the conversion of Kennicutt (1998), 
and extrapolating the flux in the U band to restframe 1500 A 
assuming a flat continuum in /„, the derived SFRs range from 
2 — 21 M Q yr -1 , with a median value of 6.0 M Q yr -1 . These 
are the SFRs uncorrected for dust, which affects the UV light 
strongly. However, as the exact dust attenuation has been deter- 
mined with the SED fitting, these can be corrected for the dust 
absorption using the same Calzetti et al. (2000) law as in the ac- 
tual SED fitting. After correction of the dust absorption, the dust 
corrected SFRs range from 5 — 64 M Q yr _1 , with a median SFR 
of 25.4 M Q yr" 1 . 



5. Discussion 

5. 1. Stellar properties of z 



1 LBGs 



The results presented here come from fitting the SEDs of high 
redshift galaxies with theoretical templates. This fitting is ham- 
pered by a number of systematic uncertainties. Firstly, on the ob- 
servational side, the inclusion of the grism spectrum in the pho- 
tometric points helps constraining the fits very well, although 
for object LBG.49138 and LBG.53929, the strong [OIII] and 
[Oil] emission lines had to be removed manually since they were 
difficult to fit accurately. One may further ask if the light from 
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these galaxies may be contaminated by e.g. AGN. We searched 
the Chandra X-ray 2 Ms catalogue (Luo et al. 2008) for coun- 
terparts of the galaxies and found only one, corresponding to 
LBG_66235. The X-ray flux of this galaxy is 1.8 x 10 42 erg s" 1 
in the 0.5 — 2 keV band, corresponding to a SFR of 355 M© yr _1 
(Ranalli et al. 2003). This SFR is much higher than that derived 
from the UV (64 M© yr _1 ) indicating that this galaxy may be a 
Seyfert galaxy. 

On the modeling side, uncertainties remain in the treatment 
of very young populations (the nebular emission) and in older 
populations (e.g. TP-AGB stars). The use of an MCMC fitting 
method reveals the blunt truth about these uncertainties. We are 
unable to fit parameters such as metallicity, SFR, and old ages. 
On the other hand, the three simple parameters of young popu- 
lation age, Ay, and stellar mass are relatively well constrained. 
This is due to the fact that the young age is mainly determined 
from the Balmer break, the Ay mainly from the UV slope, and 
the stellar mass from the overall normalisation of the SED, three 
parameters which are well constrained in the SEDs fitted here. 
Overall, the constraints on these parameters should be robust. 

The ages of the galaxies presented here are very young, and 
based on the appearance of many of them (including several 
grand spiral galaxies) it is clear that several, if not all, must have 
a longer history of star formation. Unfortunately, any older stel- 
lar populations are hidden due to the brightness of the younger 
stellar population, a bias that is true for observers of LBGs at 
both at low and high redshift. Hence, the ages presented here are 
merely a measure of the latest burst of star formation. Seen from 
the perspective of studies of LBG galaxies, the ages agree very 
well with previous results, especially those of very high redshift 
surveys from e.g. Verma et al. (2007) or Yabe et al. (2009). As 
LBGs are selected based on their bright UV colours they are 
expected to be strongly star forming galaxies in the high red- 
shift Universe. The young population ages consistently found 
in SED fitting confirms this expectation. The dust contents are 
slightly higher than the values found at high redshift; here we 
find dust extinctions of Ay » 0.5 — 2.0 mag, compared to a 
median dust extinction of Ay « 0.7 — 1.2 mag in the z ~ 2.5 
results of Shapley et al. (2001) and Papovich et al. (2001) or of 
Ay rs 0.3 mag in the z ~ 5 results of Verma et al. (2007). This 
also seems reasonable considering the potentially longer evolu- 
tion time. Masses are similar at both high and low redshift, with 
the lower redshift galaxies having a larger spread going up to 
masses of 10 11 M©. 

Overall, the properties of z ~ 1 LBGs are very similar to 
those of higher redshift LBGs. This is contrary to what has re- 
cently been seen in another large group of high redshift galaxies, 
the Lya emitters (LAEs). For LAEs, it has been reported that 
lower redshift LAEs have different properties to higher redshift 
LAEs (Finkelstein et al. 2009, Nilsson et al. 2009a, Guaita et 
al. 2010, Nilsson et al. 2010). In these publications, it has been 
shown that galaxies selected through the Lya technique at lower 
redshifts tend to be more massive and dustier than those selected 
at higher redshifts. The reason for this evolution is expected to 
be due to the fact that selecting Lya emitting galaxies will se- 
lect different objects at different times, due to the number den- 
sity evolution of objects capable of Lya emission across time. 
However, whereas Lya can be produced by many sources, the 
strength of the Lyman break of a galaxy will only depend on the 
star formation history of the galaxy and the absorption of the in- 
tergalactic medium. This could be the reason why LBGs exhibit 
similar properties across time, as opposed to LAEs. Hence, the 
selection methods finding LBGs and LAEs are in some ways or- 
thogonal, in that the LBG technique will find the same galaxies 



at all redshifts, and the LAE method will find a representation of 
all Lya emitting objects in each redshift slice it surveys. 

5.2. The importance of the grism photometry 

To test the importance of the grism spectral points in the fits, 
we also ran the SED fits for various combinations of photo- 
metric points included. For this test, LBGs_# 59257, 63984, 
66235 and 69919 were selected as test cases as they represent 
a range in x 1 , morphology and number of grism points. The 
SEDs were divided into four parts; the UV/optical part ("O"), the 
grism points ("G"), the near-infrared points ("N") and the single 
Spitzer/IRAC Chi point ("S"). New SED runs were then made 
with realistic permutations of these photometric parts. The re- 
sults for one object, LBG_69919, are presented in Fig. [2] For the 
analysis of these results, the values of the fit including all bands 
(OGNS) are assumed to be the most correct, standard value. In 
Fig. |2] both the accuracy of the fit, in comparison to the OGNS 
result, and the significance of the result are shown. Points within 
the solid line cone all have parameter results that agree with the 
OGNS result within la. Fits without the grism photometry in- 
cluded are all shown with green symbols. 

Several conclusions can be drawn from this Figure and the 
results of the fits of the other three objects. Firstly, most of all fits 
including the grism points are within la of the OGNS values, 
whereas the fits without the grism points often lie far from the 
standard values. Thus, the grism points are necessary to get the 
best possible result. Another conclusion, and also part of the rea- 
son why the fits with grism points are always within la, is that 
the error bars are typically larger on the fits with grism points, i.e. 
a lower significance of the parameter. This can be slightly sur- 
prising, considering that the x 2 r are consistent throughout all fits 
(e.g. it varies between 0.8 — 3.1 for LBG_69919). One possible 
reason for this is that the stellar templates used are not accurate 
enough in predicting the Balmer break in the galaxies. With up 
to 68 photometric points determining the Balmer break, albeit 
with large error bars, small uncertainties in the spectral popula- 
tion modeling could potentially cause large uncertainties in the 
fitted parameters. 

Further, fits including the Spitzer/IRAC Chi point are typ- 
ically better fit than those without. An important conclusion is 
thus that it is always important to sample also the restframe near- 
infrared when fitting stellar populations. Finally, in Fig. [2] it is 
seen that fitting the SEDs consisting of only the grism points re- 
veals parameter results that are within la envelope of the OGNS 
results. The significance of the results in this case are not very 
high, but the parameters are constrained enough to give inter- 
esting results, also in studies where other multi-wavelength data 
is not available. This conclusions seems to break down though, 
as fewer grism data-points are available, as in the case of e.g. 
LBG.59257 or LBG.66235 . 

Note though, that the conclusions presented here are valid 
for galaxies at z ~ 1 where the grism points sample the Balmer 
and 4000 A breaks of the galaxies. It is unlikely that fitting 
only grism points would give good results for other redshifts, 
although it cannot be excluded from the present data. 

5.3. Conclusion 

To conclude, we have demonstrated that including grism spec- 
trophotometric points in SED fitting of z ~ 1 LBGs have im- 
proved the fits of these galaxies greatly. The z ~ 1 LBGs all 
have very young ages, possibly with one or more underlying 
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older populations outshone from the brighter younger popula- 
tion. In order to be able to constrain older populations better, 
more accurate photometry in the restframe near-infrared bands 
(observed Spitzer bands) would be necessary, which is challeng- 
ing for galaxies at z ~ 1 and even more difficult for higher red- 
shift galaxies with existing observatories. The LBGs have mod- 
erate to high dust extinction and the masses are similar to those 
of higher redshift LBGs. Considering the similarity of properties 
of LBGs at higher and lower redshifts, it appears that selecting 
high redshift galaxies according to their UV properties selects 
the same type of galaxy, irrespective of redshift, as opposed to 
the selection via e.g. Lya emission. 

References 

Balestra, I., Mainieri, V., Popesso, P., et al., 2010, A&A, 512, A12 
Balogh, M.L., Moms, S.L., Yee, H.K.C., Carlberg, R.G., & Ellingson, E., 1999, 
ApJ, 527, 54 

Bolzonella, M., Miralles, J.-M., Pello, R., 2000, A&A 363, 476 
Brazual G.A., Chariot S., 2003, MNRAS, 344, 1000 

Bunker, A.J., Stanway, E.R., Ellis, R.S., & McMahon, R.G., 2004, MNRAS, 
355, 374 

Burgarella, D., Le Floc'h, E., Takeuchi, T.T., et al., 2007, MNRAS, 380, 986 
Calzetti, D., Armus, L., Bohlin, R.C., et al., 2000, ApJ, 533, 682 
Daddi, E., Cimatti, A., Renzini, A., et al., 2004, ApJ, 617, 746 
Finkelstein, S.L., Cohen, S.H., Malhotra, S., & Rhoads, J.E., 2009, ApJ, 700, 
276 

Franx, M., Labbe, I., Rudnick, G., et al., 2003, ApJ, 587, L79 
Fynbo, J.P.U., M0ller, P., Thomsen, B., et al., 2002, A&A, 388, 425 
Gawiser, E., Van Dokkum, P.G., Herrera, D., et al., 2006, ApJS, 162, 1 
Giavalisco, M., Ferguson, H.C., Koekemoer, A.M., et al., 2004, ApJ, 600, L93 
Guaita, L., Gawiser, E., Padilla, N., et al., 2010, ApJ, 714, 255 
Gunn, J.E., & Stryker, L., ApJS, 52, 121 
Hamilton, D., 1985, ApJ, 297, 371 

Kauffmann, G., Heckman, T.M., White, S.D.M., et al., 2003, MNRAS, 341, 33 
Kennicutt, R.C., 1998, ARAA, 36, 189 

Kriek, M., van Dokkum, P.G, Franx, M., et al., 2006, ApJ, 645, 44 

Le Fevre, O., Vettolani, G., Garilli, B., et al., 2005, A&A, 439, 845 

Leitherer, C, Schaerer, D., Goldader, J.D., et al., 1999, ApJS, 123, 3 

Luo, B., Bauer, F.E., Brandt, W.N., et al., 2008, ApJS, 179, 19 

Mignoli, M., Cimatti, A., Zamorani, G., et al., 2005, A&A, 437, 883 

M0ller, P., & Wan-en, S.J. 1993, A&A 270, 43 

Nilsson, K.K., M0ller, P., Moller, O., et al., 2007, A&A, 471, 71 

Nilsson, K.K., Tapken, C, M0ller, P., et al., 2009a, A&A, 498, 13 

Nilsson, K.K., Moller-Nilsson, O., M0ller, P., Fynbo, J.P.U., & Shapley, A.E., 

2009b, MNRAS, 400, 235 
Nilsson , K.K., Ostlin , G., M0ller, P., et al., 2010, submitted to A&A, 

arXiv: 1009.0007 

Nonino, M., Dickinson, M., Rosati, P., et al., 2009, ApJS, 183, 244 
Ouchi, M., Shimasaku, K., Okamura, S. et al., 2004, ApJ, 611, 660 
Overzier, R.A., Heckman, T.M., Tremonti, C, et al., 2009, ApJ, 706, 203 
Papovich, C, Dickinson, M., & Ferguson, H.C., 2001, ApJ, 559, 620 
Pettini, M., Shapley, A.E., Steidel, C.C., et al., 2001, ApJ, 554 981 
Pirzkal, N., Xu, C, Malhotra, S., et al., 2004, ApJ, 154, 501 
Ranalli, P., Comastri, A., & Setti, G, 2003, A&A, 399, 39 
Ravikumar, CD., Puech, M., Flores, H., et al., 2007, A&A, 465, 1099 
Retzlaff, J., Rosati, P., Dickinson, M., et al., 2010, A&A, 511, A50 
Shapley, A.E., Steidel, C.C., Adelberger, A.E., et al., 2001, ApJ, 562, 95 
Shapley, A.E., Steidel, C.C., Pettini, M., & Adelberger, K.L., 2003, ApJ, 588, 65 
Steidel, C.C., Giavalisco, M., Pettini, M., Dickinson, M., & Adelberger, K.L., 

1996, ApJ, 462, L17 
Steidel, C.C., Adelberger, K.L., Giavalisco, M., Dickinson, M., & Pettini, M., 

1999, ApJ, 519, 1 

Straughn, A.N., Pirzkal, N., Meurer, G. R., et al., 2009, AJ, 138, 1022 
Vanzella, E., Cristiani, S., Dickinson, M., et al., 2008, A&A, 478, 83 
Venemans, B.P., Rottgering, H.J.A,; Miley, G.K., et al., 2007, A&A, 461, 823 
Verma, A., Lehnert, M.D., Forster-Schreiber, N.M., Bremer, M.N., & Douglas, 

L., 2007, MNRAS, 377, 1024 
Xu, C, Pirzkal, N., Malhotra, S., et al., 2007, AJ, 134, 169 
Yabe, K„ Ohta, K„ Iwata, I., et al., 2009, ApJ, 693, 507 



o 

(0 



a 

60 




L 
0.0 



12 



0.5 1.0 1.5 2.0 

Age of fit compared to standard 



2.5 



10 



m 

















OGNS • 










G • 










GN 










GS • 






* 1 




GNS 










OG 










OGN 










OGS • 






1 













\ * 


ON 






1 




OS 






I 




ONS 








• \ 








• \ 








• \^ 













0.0 0.5 1.0 1.5 

A v of fit compared to standard 



2.0 



20 



o 15 



10 



a 

at) 
m 

























OGNS • 






: I 






O 


G • 






• 1 




i 




GN 










I ; 




GS • 














GNS 






: t 








OG 














OGN 














OGS • 





























ON 














OS 














ONS 






/ 
























m 


• 

















0.0 0.5 1.0 1.5 

M. of fit compared to standard 



2.0 



Fig. 2. Test of the importance of various photometric bands for 
the fitting. The three plots are shown for the parameters age, 
Ay and stellar mass, from top to bottom. The x-axes show the 
value for each fit normalised to the value obtained with all fil- 
ter passbands used (OGNS, black points). The y-axes show the 
value in each fit divided by the error bar on that value, indicat- 
ing how well determined the value is. The various passbands 
included are UV/optical ("O"), grism ("G"), near-infrared ("N") 
and Spitzer/Chl ("S"). The solid, dashed, and dotted line indi- 
cates the area within which points are within 1, 2, and 3<r of the 
OGNS value. 
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Fig. 3. Typical dependencies between parameters. From the 1 SSP fits to LBG_46190. Contours show la, 2a and 3a confidence 
levels. The green point with error bars are the best fit ID values, as given in Table[2] The metallicity is completely unconstrained, 
whereas the parameters Ay, stellar mass and age are relatively well constrained. There is a slight degeneracy between age and dust, 
with older ages requiring less dust. 



